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We analyze an approach aiming at determining statistical properties of spectra of time-periodic 
quantum chaotic system based on the parameter dynamics of their quasienergies. In particular we 
show that application of the methods of statistical physics, proposed previously in the literature, 
taking into account appropriate integrals of motion of the parametric dynamics is fully justified, 
even if the used integrals of motion do not determine the invariant manifold in a unique way. The 
indetermination of the manifold is removed by applying Dirac's theory of constrained Hamiltonian 
systems and imposing appropriate primary, first-class constraints and a gauge transformation gen- 
erated by them in the standard way. The obtained results close the gap in the whole reasoning 
aiming at understanding statistical properties of spectra in terms of parametric dynamics. 
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I. INTRODUCTION 



One of the most characteristic features of quantum systems which exhibit chaotic behaviour in the classical limit 
. is an affinity of their spectral properties to random matrices. The famous Bohigas-Giannoni-Schmidt conjecture 
states that the statistics of distances between neighbouring energy levels of a quantum system with chaotic classical 
limit is well described by the one derived from the Random Matrix Theory (RMT) Q ■ 

A vast numerical and experimental evidence 0, 0] in favour of this hypothesis was collected during last twenty 
years. There are also convincing theoretical arguments supporting it lajfj- In the present paper we would like to 
reconsider one of the first theoretical approaches initiated by Pechukas and further developed by Yukawa 
The original idea consisted in deriving differential equations describing parametric level dynamics i.e. the evolution 
of eigenvalues, when the parameter controlling the amount of chaos in the system changes, and applying the rules 
of classical equilibrium statistical mechanics to the flow described by the derived differential equations, treating the 
^ parameter as a fictitious time in which the "evolution" takes place. As observed by Yukawa the resulting dynamical 

■ system was Hamiltonian, hence applying rules of equilibrium statistical mechanics was straightforward; the equilibrium 
OO ' distribution should be given as the Boltzmann one, 

o ' 

VO ■ p = Afexp(-PH), (1) 

O . 

where Ti. is the Hamilton function of the system, - a fictitious temperature (to be determined in some way from 
the initial data), and M - an appropriate normalization constant. The dynamical variables of the model, apart from 
the energy eigenvalues, involved also other ones. Integration over them over the phase space led to the equilibrium 
distribution of energy levels. As shown by Pechukas and Yukawa the resulting distribution coincides with those 
provided by RMT for the ensemble of real symmetric matrices with identically and independently distributed elements 

■ (forming the so-called Gaussian Orthogonal Ensemble). 

As appealing and straightforward as the above outlined approach might be, one should not overlook some funda- 
mental obstacles appearing when attempting to formulate it in a more rigourous way. Let us summarize briefly the 
most disturbing of them. 

Both Pechukas and Yukawa started with the quantum Hamiltonian of the form 

H = H + XV, (2) 

where, according to their original interpretation, a time-independent N x N Hermitian matrix Hq represented a 
quantum system enjoying integrable classical limit, whereas V was an integrability-breaking part making the whole 
system classically fully chaotic when A attained appropriately large values. Thus the ambitious program was designed 
to actually investigate the transition between spectra of integrable (A = 0) and nonintegrable (A - large) cases. Thus, 
treating A as a fictitious time, one faces a problem belonging to non-equilibrium rather than equilibrium statistical 
mechanics, and usefulness of tools of the latter could be doubtful. Moreover, as it is clear from (J2J, as A grows 
the motion is unbounded, in particular the eigenvalues of H grow indefinitely and, without an additional scaling, 
no 'equilibrium' distribution of eigenvalues of H is attained (although it might be that the statistics of distances 
measured in units of the mean distance approaches some 'equilibrium'). There is a way of curing the situation - one 
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should turn to dynamics of rescaled energy levels, what in fact consist in changing the A-dependence of H (see Q). In 
the following we take another route - we would like to investigate the spectra of unitary evolution propagators instead 
of Hermitian Hamiltonians. Propagators being unitary have their spectra confined to the unit circle independently 
of their specific parameter dependence. The idea goes back to Dyson [T^j, although in slightly different context - he 
realized that it is more convenient to define probabilistic measures on ensembles of unitary rather than Hermitian 
matrices, what eventually led to the definitions of the circular ensembles of RMT. 

In our case we do not pretend to mimic the same way. Instead we propose to start with the propagator for a 
particularly simple time-dependent quantum system in which the integrability breaking part in J2J has a form of 
periodic instantaneous kicks, so the whole Hamiltonian reads now: 

oo 

H(X) = H + XV S(t-nT), (3) 



where V is some constant, Hermitian matrix and T the period of the perturbation. 

The unitary evolution operator of the system (the propagator) is given as the solution of: 

HF 

ih— = H(X)F, F(0)=I. (4) 

Since the time dependence is periodic the whole information about the evolution of the system is encoded in F 
evaluated at time T, i.e. the propagator transporting the system in time over one period of the perturbation. Since 
we will be concerned with properties of only this particular period-one propagator we will use F to denote it without 
risking confusion with the general time-dependent one. In case it is needed we will write F(X) to remind its dependence 
of the perturbation parameter A inherited from the original Hamiltonian • In the following we put h = 1 and T = 1. 
In the case of a kicked system, the one-period evolution operator F takes a particulary simple form 

F(X) = exp {-iXV) F , F := exp(-iH ). (5) 

Of our interest will be the eigenphases ip n {X) (i.e. the phases of the eigenvalues), sometimes called also quasienergies 
ofF(A), 

F(A)|0„(A)) = exp(^„(A))|^(A)), (6) 

where \(j) n (X)} are eigenvectors of F(X). 

Symmetries of the system in question (in particular with respect to the reversal of time) enforce additional symmetry 
conditions on H and F. Without any other symmetries present, H is a general hermitian and F a general unitary 
N x N matrix (we will be mostly concerned with this case), whereas eg. for integer-spin time-reversal-symmetric 
systems H is real symmetric H* = H = H T and we can assume F = F T (i.e. F symmetric). 

As in the autonomous case (0) considered by Pechukas and Yukawa one can derive a closed set of differential 
equations, describing a Hamiltonian motion in a multi-dimensional phase-space with A as a fictitious time ( HGH ) 
(see next section). Since eigenvalues of an unitary matrix lie on the unit circle, their phases are restricted to an 
interval [0,27r[, and in contrast to the autonomous case (0), the motion of the relevant variables is now bounded. 
This observation eases to some extend the above mentioned problems with the Pechukas- Yukawa approach, although 
does not warrant any equilibration as A grows. The problem of attaining the equilibrium was thoroughly discussed in 
|13| where it was shown that the time of effective equilibration goes to zero when approaching classical limit of the 
systems in question, hence the equilibrium statistical properties of the spectra should be detectable for large quantum 
numbers. We will recapitulate this discussion in the following sections. 

When passing to the distribution of eigenphases from the equilibrium distribution involving all dynamical 
variables one should choose an appropriate measure for integrating out the irrelevant variables. The natural measure 
is inherited from the symplectic structure underlying the parametric dynamics of the eigenphases and was discussed 
in 0. 

Finally let us introduce briefly the main objective of the present paper aiming at closing the final gap in vindication of 
the parameter dynamics approach to statistical properties of quasi-energies outlined above. The postulated Boltzmann 
distribution can be validated only if there are no other constants of the motion apart from the Hamilton function itself, 
or in other words, when the motion is ergodic on the whole constant-energy surface. In the case of Pechukas- Yukawa 
parametric dynamics, as well as in the case of parametric motion of the kicked periodic systems © it is not the 
case. The dynamical systems governing the parametric motion of eigenvalues or eigenphases are so called generalized 
Calogero-Moser or Sutherland-Moser systems . They posses many additional integrals of motion and in fact the 
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motion is ergodic on a much smaller invariant manifold [L6(. The simplest way to include the influence of additional 
integrals of motion consists of using in place of canonical ensemble measure Q its grand-canonical generalization 

p oc expj - ^Pn 1 ^} ; ( 7 ) 

nailing down the invariant manifold on which the motion takes place by fixing constants of the motion / M in the 
ensemble mean with the help of Lagrange parameters (3^; one of these I n should be the Hamilton function 7i. Using 
the microcanonical ensemble p oc Y\ ${In fixing the values of the constants of motion I n to their initial values I n 
exactly, rather in the ensemble mean as Q does, would be even more appropriate, but technically more complicated 
- see [3j , Chapter 6 for a discussion of the problem. 

Integration of p over all dynamical variables except the eigenphases yields the desired distribution P(tpi, ■ ■ ■ , <Pn)- 
Such a program was performed in [171 Il8| where it was shown that inclusion of additional known integrals of motion 
leads to corrections of the order X/N in comparison with the predictions of RMT. This result is highly satisfactory, 
since one expects convergence to RMT in the limit when the dimension of the matrix N tends to infinity (which, 
in many models corresponds to classical limit of the quantum system). The only remaining problem is whether the 
integrals of motion taken into account in [Til Hsj are all, which are needed to fix (in the ensemble mean) the invariant 
manifold on which the motion is ergodic, or in other words, what is the minimal set of independent integrals of motion 
determining the invariant manifold (see also Q for the formulation of the problem). The independence of the integrals 
used in the above mentioned papers was investigated in |19|. In the present paper we close the last gap by showing 
that they determine the invariant manifold to the extent which is appropriate to infer the distribution of eigenphases. 

In Section [n] we briefly derive the dynamical equations for the parametric motion of eigenphases, whereas in 
Section IIIII we discuss the statistical mechanics of the one-dimensional gas of eigenphases. Section IIVI is devoted 
to additional integrals of motion: we give their form appropriate for further applications and briefly discuss their 
independence. The final Section Ivl offers the ultimate solution of the problem in frames of the theory of constrained 
Hamiltonian systems. 



II. PARAMETRIC LEVEL DYNAMICS 



Let W(X) diagonalizes F (A), 

W(X)F(X)W- 1 (X) = e -**( A ) =diag(e-^ I W,...,e-^W), (8) 

$(A)=diag^i(A),...,^jv(A)). (9) 

In the following we shall skip exhibiting the explicit A-dependence when possible. Since F is unitary, so is W, 
W~ x = W'. In the case of a general unitary matrix, the diagonalizing matrix W is not unique even after ordering 
the eigenphases: we can always left-multiply it by a diagonal unitary matrix without altering the result (JSJ). 
Let's define following auxiliary matrices 

v := WVW' 1 =v\ (10) 
I := ie l * he" 8 *] = -Jt. (11) 

^From Ijllfl and the diagonal character of e _I * we have 



Inn = 0, V mn = - -r^— , U ^ TO. (12) 



Differentiating ijSJ over A we arrive at 



dA 

dv 

dX 
dl 

dX 



where 



i(a-e- l *ae l *)+v, (13) 
[a,v], (14) 
[a,l], (15) 



« = ^-U (16) 



dX 
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For reasons which soon will be clear in the following we shall use the notation: 

q n := tp n , p n := V nn , (17) 
The diagonal part of the matrix equation 1)13)1 reads 

-jY = v nn =p n , (18) 
whereas its off-diagonal part gives the off diagonal elements of a in terms of v mn 



what, upon (|T2|) . gives: 

I, 



, to ^ n, (19) 



Qmn = —. !""o „ , m^n. (20) 
4 sin z - 



By an appropriate choice of the diagonalizing matrix W - see remark, below (JHJ), we can choose a nn = 0. 

Eliminating with the help (|12f> and 1)19)1 v mn and a mn in favour of l mn we obtain from the diagonal part of (|14fl 



= l nkhnV'(qk - q n ), (21) 



and, form 1)15)1 for n ^ to, 



where 



M = - >, ^mfc'fen (V(g„ - - V(qk - q m )) , (22) 



4 sm z | 



and ' in l|21|l denotes the derivative with respect to the argument. 

Equations (|18|) . I|21|) . and (1221) are Hamiltonian with A treated as a fictitious time 



| = {«,!»}, ^T = Wl>J> ^f = W m J, (24) 



and with the Hamilton function 

N , iV 



^ = \ ^jPI + \ l mJn m V(q n - g m ) = ^ Tr-y2 ' ( 25 ) 
n— 1 71,771=1 

if we define the following Poisson brackets among the dynamical variables q n ,Pni and lmn 

{Pm, q-n} = <W) {Pm.Pn} = {<7m, <7n} = 0, (26) 

IP™, ^fcn} = {9m, Jfcn} = 0. (28) 



It might be appropriate to mention that the system of dynamical equations for the parametric motion of eigenvalues 
of a Hermitian matrix considered by Yukawa [gj has the same form as 1)181) , 1)21)1 , and 1)22(1 but with a different form of 
the potential, V(q) = —1/q 2 , the 'spatial' coordinates corresponding to eigenvalues, and slightly different definitions 
of v and I. For a unified treatment of parametric motion in autonomous and kicked cases see |l ll l20j . 

The system of equations l(18|) . 1)21 [I. and 1)22)1 can be treated as describing dynamics (in the fictitious time A) of 
a one-dimensional gas of particles on the unit circle interacting mutually via the potential 1)23)1 . but with evolving 
'coupling strengths' l mn becoming thus additional dynamical variables. 
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III. STATISTICAL MECHANICS OF THE GAS OF EIGENPHASES 



As the phase-space trajectory of the fictitious gas evolves in time A, the original matrix F(X) changes within a one- 
parameter family. It is that family rather than a single dynamical system (which has a fixed value of A) which exhibits 
random-matrix type spectral fluctuations. Indeed, if the motion is ergodic (we shall show below, that it is indeed, 
although not on the whole energy surface) implies that A averages of spectral characteristics like the distribution 
of spacings between adjacent quasienergy levels equal ensemble averages. Of importance is thus the minimal time 
interval AA needed for time and ensemble averages to become practically equal. Obviously AA — > oo is sufficient, 
but not necessary provision. During the evolution 'particles of the gas' (i.e. in fact, the eigenphases) undergo mutual 
collisions. Observe that since the potential in l|25(l is repulsive, they usually do not cross (i.e. do not exchange 
positions). Such a real crossing of two eigenphases would demand vanishing the respective value of l mn . Instead 
what is usually observed in the region of parameter A corresponding to classically chaotic behaviour, are so called 
avoided crossings when two neighbouring quasienergies approach a minimal nonzero distance when A changes. In 
fact, as numerical experiments show, for systems which are classically chaotic such avoided crossings are abundant 
[l3|. Due to collisions the gas reaches state in which the motion of particles consists of fluctuations in the vicinity of 
equilibrium. If Hq is integrable and chaos develops gradually after switching the perturbation V and increasing the 
coupling strength A, the spectrum of F(X) equilibrate to RMT predictions only after certain 'relaxation time' when 
the phase space regions of regular motion have shrunk to relatively negligible weight. If, on the other hand, as shown 
m lig, Ho and V are both non-integrable the initial state of the fictitious gas is already close to equilibrium, and then 
the window AA in question need not be much larger than the collision time A co ii of the gas, i.e. the mean distance of 
avoided crossings for a pair of neighbouring levels. It was shown in |13| that the time elapsing between consecutive 
collisions scales as N~ v , v > 0. A A-average for the over a window AA oc A co n thus involves a family of operators 
F(X) which all yield identical classical dynamics in the limit N~ u — > 0. 

The most straightforward application of statistical mechanics is to employ the canonical ensemble for the distribution 
of the dynamical variables (q,p, I) 

p(q,p,l)<xe X p{-W(q,p,l)). (29) 

A straightforward integration over Gauss-distributed p and I gives precisely the eigenphase density of random-matrix 
theory 0, i.e. 

P(q u ...,q N ) = [dVp^-Wle-f™* J[ | e -*» - e -^| (30) 

J m<n 

As explained in the introduction the reasoning would be reasonable, if there were no other integrals of motion beside 
TL itself. In the case other integrals of motion 1^ exist, the appropriate ensemble to use is the generalized canonical 
ensemble (JTJl. 



IV. INTEGRALS OF MOTION 



The equations H18|) . I21|) . and ll'L'li are clearly integrable (they can be solved simply by diagonalizing F at given 
A and calculating appropriate matrix elements), so one should expect that there are much more integrals of motion 
than the Hamilton function l|25|) itself. Indeed from i|13fl and 11141 we see that the quantities 

h imi ...k nmn = r Tr(v k H m ^--v k -r"), (31) 

are indeed constants of motion, i.e. 

dX 



and should be taken into account when constructing the generalized canonical ensemble (J7J). It can be shown |3l ll7lll8| . 
that such an ensemble yields the distribution of level spacings as well as low-order correlation functions of the level 
density in common with random-matrix theory, to within corrections of order 1/N. The only problem is whether all 
integrals nailing down an invariant manifold are of the form (I31|l . 

It is possible to show that only N 2 — N of such integrals are independent. Indeed, let us briefly recall the reasoning 
presented in |l9j and instead of Ik 1 m 1 ...k Tl m Tl consider N 2 quantities 

C km := Tr (e**t> fc e- i *« m ) , 0<k,m<N-l. (33) 



6 



^From the definition of I (|llfl 



= v — il, 



(34) 



i.e. 



e iZ> v k e -i<S> v m 



(35) 



(v-ilfv m , 
hence Ckm are linear combinations of l|31|) . 

1. Not all Ckm are independent. Indeed, trivially Cko = Cot = TrV k . 

2. There are no more independent integrals of this type (i.e. with k > N or m > N). Indeed from the Cayley- 
Hamilton theorem applied to the matrix v, a Ckm with k or m larger than N — 1 can be expressed a linear 
combination of the basic ones l)33|) . 



Using JHJ), and (|10[) we can rewrite Ckm as 

Ckm 



= Tr (WFWv k WFW*v m ) 
= Tr (F^(WU k W)F(WU m W)) 
= Tr (F^V k FV m ) 

= Tr \e- iXV U V k Ule iXV V™ 
= Tr ( U a V k UlV r 



(36) 



In the basis in which V is diagonal i.e. Vij = Vrfij this reduces to 

N 

Ckm = £ (V,) fc (V p r|(^ ) OT | 2 , 



< k,m < N - 1. 



(37) 



Now l)33|l can be treated as a system of TV 2 linear equations for TV 2 unknown quantities |(f/o)pg| 2 j 1 < Pi Q < AT, to 
be expressed in terms of A^ 2 quantities Cfc m , < fe,m < JV — 1. For a generic V i.e. when all its eigenvalues Vi are 
different the system can be solved, since the determinant of the coefficient matrix is a power of the Vandermondc 
determinant Dy constructed from Vi, 



D 



v •= 



1 Vi 

1 V2 

1 Vat 



V, 



N-l 
'N—l 



V, 



N-l 
N 



(38) 



K 3 



Thus all Ckm can be expressed as linear combinations of \(Uo) pq \ 2 . The iV 2 real numbers Uij := \(Uo) pq \ fulfill 2N 
relations stemming from the normalization of rows: 



and columns 



JV 
3=1 



N 



i, i = i, 



(39) 



(40) 



of the unitary matrix U$. The number of independent relations is equal to 2N — 1 and is less by one than the 
total number of equations in l|39|l and 14U|) since summing all equations in l|39|l over i gives the same as summing all 
equations in |@DJ over j, namely N = TtUqUq. Finally thus all \(U ) pq \ 2 involve N 2 - (2JV-1) = (N-l) 2 independent 
parameters and the number of independent Ckm which can be used to determine them must be at least equal. Since 
Cko and C fc do not depend at all on Uo, they can not be used to determine \(Uo) pq \ 2 . To do this we are left only 
with Ckm, 1 < k,m < N which are exactly (N — l) 2 in number. The Uq - independent integrals Cko, < k < N, on 
the other hand, can be used to determine N independent parameters of V (traces of its powers, or what is equivalent, 



7 



its eigenvalues), so they are also independent. Since, trivially, Coo = N, we are left with the independent integrals of 
motion of the form: 

C k :=C k0 = TrV k , k = 1, 2, . . . , N - 1, (41) 
which are TV — 1 in number, and (N — l) 2 integrals 

C km = Tr (u Q V k UlV m ) , k, m = 1, 2, . . . , JV - 1, (42) 

i.e. all together N 2 — N independent integrals of motion. 

Let us now count how many variables we have in l)18[l. 121|l . and Q22JI. The variables q n are real as eigenphases 
of the unitary matrix F and there are N of them. Also p n as diagonal elements of a Hermitian matrix v are real, 
there are N of them as well. Since I is antihcrmitian and off-diagonal there are (N 2 — N)/2 matrix elements l mn , but 
since they are, in the case of a general unitary matrix, complex, we should count separately their real and imaginary 
parts. Finally thus we have N + N + (N 2 — N) = N 2 + N real variables. Comparing this with the number of found 
integrals of motion (./V 2 — N) we are tempted to think that, in a generic case, invariant manifolds are of the dimension 
(N 2 +N)- (N 2 — N) = 2N. 

On the other hand, in the coordinate frame in which V is diagonal, = ViSij, the motion described by F(X) = 
exp(— iV)Uo involves only N independent frequencies Vfc, and takes place on an iV-dimensional torus (and is ergodic 
on it i a generic case when the eigenvalues of V are not rationally dependent). It seems thus, we are still missing N 
independent integrals of motion. 



V. POISSON STRUCTURE AND CONSTRAINTS 



Before identifying missing integrals and determining their influence (or lack of) on the distribution of eigenphases, 
let's consider more carefully the proposed Hamiltonian formulation. First, observe that the definition of the manifold 
on which the level dynamics takes place as parameterized by the coordinates q n ,p n and l mn and equipped with the 
Poisson structure l|26(l - (|28|l is slightly flawed. From the definition (jl 1|> of I we have l nn — 0, but this is inconsistent 
with the Jacobi identity which must be fulfilled by l|28|l : 

Indeed, substituting to the above m = k and n = i and using =1^ = whenever they appear on intermediate 
steps, we arrive at an erroneous result S pq lk q + 5i q l P i — Skqlpk — Spihq- instead of zero. Thus we are not allowed to put 
Inn = from the very beginning as equations defining our manifold. Instead, if we want to keep the Poisson brackets 
(|28H we should change the definition fTJ to 

l = ie l * [v^e-^] +iL, (44) 

where L is an arbitrary, real, diagonal matrix, i.e. we introduced N additional dynamical variables. To understand 
their meaning let us return to the derivation of the dynamical equations by diagonalizing matrix W IjSjl. but this 
time we do not impose additional conditions on W, i.e. we do not assume that the diagonal matrix elements a nn of 
a = dW/dX ■ W~ x vanish. Instead we allow them to be arbitrary functions of A. It should be clear (and indeed we 
will show that it is the case), nothing really depends on the choice of a nn , since nothing concerning the eigenvalues 
should depend on the choice of the diagonalizing matrix. 

The resulting equations of motion are derived in the same way as previous ones JTSJl, (|21|l . and (|22|l . In fact only 
the third of them is altered and reads now, 

= - Imkhn (V(g„ - <?fc) - V{q k - q rn )) 

+ -Qn)- (45) 

Equations (|18fl . I|21|l . I|45|) are again Hamiltonian with the same Poisson structure (|26l) - (|28|) . but with a new Hamilton 
function 

N N N 

n = 2^2 p2n+ 2 ^2 l ™nl nm V(q n - q m ) + ^ a n l ii > ( 46 ) 

ri—l n.m— 1 j 
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depending on N arbitrary (in general 'time-', i.e. A-dependent) functions a nn . The quantities C mn (|33|l are again 
integrals of motion. In addition, we easily calculate that 

{H, Inn} = 0, (47) 

so Inn are also constants of motion. In fact, as it is clear from the previous considerations, nothing concerning the 
eigenphases depends on actual values of l nn . We can thus impose constraints, eg. 

l nn = 0, n = l,...,N. (48) 

At this point it is instructive and in fact very natural to describe the encountered situation from the point of view 
of Dirac's theory of constrained Hamiltonian systems |2l|. The conditions (|48|l are so called primary constraints (i.e. 
they are not obtained form the equations of motion) and can be imposed only after evaluating all Poisson brackets to 
avoid the problems with the Jacobi identity mentioned at the beginning of the present section. Further, one calculates 
easily that 

{Immilnn} 0. (49) 

Together with (|47Jl it means that the consistency condition 

N 

lnim\ ~t~ ^ Qnn\lnni ^mm\ 0, (^0) 
n=l 

are identically fulfilled and no other constraints, neither primary nor secondary, are produced, nor additional conditions 
are imposed on the functions a„„(A). 

Due to (|49|l l nn are automatically first-class constraints (recall that according to Dirac terminology a quantity is 
of first-class if its Poisson brackets with all constraints vanish, see [2l|, p. 18). The new Hamilton function H4tj|) 
involves as many arbitrary functions (in our case these are functions a nra (A)), as there are independent primary first- 
class constraints. On the other hand, first-class primary constraints (|4*%|) may be always used to produce a gauge 
transformation generated by 

N 



G(A) = 5>„(A)Z„ n , (51) 



Itf-e^V-^W., (52) 

with arbitrary A-dependent k = 1, . . . , N. 

We expect that an initial physical state determined by initial values of the phase-space variables (q n ,Pn,lmn) 
determines also its all future physical states. Since the Hamilton function (|46J) depends on N arbitrary functions, 
the same may happen to the values of (q n ,p n , Imn) at latter times. But the only freedom is now given by the gauge 
transformation i|52|) connecting the variables describing the same physical state of the system for different choices of 
the gauge. Hence particular physical properties of the state (eg. statistical properties of the distribution of positions, 
ie., in our case, eigenphases) should be gauge- independent, and in fact they are, since the gauge transformation does 
not influence the relevant variables q n . 

To be even more concrete in explaining the role of the gauge transformation for the present problem let us observe 
that by assuming l mm — we recovered the previous count of the number of variables vs. dimension of the invariant 
manifold, since the number of variables was first increased by N by introducing the diagonal elements of I and then 
decreased by the same number by imposing constraints equating them to zero. To fix the (still) remaining N degrees 
of freedom let us observe that the gauge transformation (|52(l does not change the integrals of motion (in particular 
the Hamilton function itself) after reducing to the manifold determined by the constraints 1)48(1 . retaining also the 
equations of motion in their original form. The transformation is intimately related to the freedom of choice of the 
diagonalizing matrix W in terms of a, it leads to 

a^i^- + e w ae-' w , 6 := diag(6 u . . . , 6 2 ). (53) 
dX 

With the help of (|52|l we can fix in an arbitrary way N (more precisely N — 1, but one additional is determined by 
the choice of initial point on the unit circle) phases of the variables l nm - Let us summarize 
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• number of variables: N var = N 2 + 2N (the old ones plus the (imaginary parts of) diagonal elements of I 

• number of independent integrals C mn : Ni nt = N 2 — N 

• number of constraints l nn = 0: N c = N 

• number of phases fixed by choosing a gauge N g = N, 

hence N var — (Ni nt + N c + N g ) = N = dimension of the invariant manifold. 

Now it is clear that integrals of motion C mn (I33|l are the only quantities which should be taken into account when 
determining the equilibrium distribution. Indeed, as already mentioned the constraints and the gauge, involving only 
Imm do not influence eigenphases, what is a direct consequence of the independence of the eigenvalues on the choice 
of the diagonalizing matrix. Moreover, our choice l nn — reduces the Hamilton function l|4rj|) to originally considered 
one (|25[1 and the whole reasoning which led, after integration out of p and I variables and neglecting corrections of 
order 1/N, to random matrix results for the eigenphases, is fully vindicated. 
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